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1 Introduction 


The  NIST  Sparse  BLAS  (Basic  Linear  Algebra  Subprogram)  library  provides  computational  kernels  for 
fundamental  sparse  matrix  operations: 


® sparse  matrix  products, 

C <-  q A B + 0C 
C <r-  a At  B + PC 

• solution  of  triangular  systems, 

C <-  aDL  A~l  DRB  + 0C 
C <-  aDh  A~t  DrB  + 0C 


where  A is  sparse  matrix,  B and  C are  dense  matrices /vectors,  and  Dl  and  DR  are  diagonal  matrices. 
This  version  of  the  NIST  Sparse  BLAS  supports  the  following  sparse  formats:  compressed  sparse  row 
(CSR),  compressed  sparse  column  (CSC),  coordinate  (COO),  block  sparse  row  (BSR),  block  sparse 
column  (BSC),  block  coordinate  (BCO)  and  variable  block  row  (VBR).  Symmetric  and  skew-symmetric 
versions  are  also  supported. 

The  routines  are  written  in  ANSI  C and  are  callable  from  Fortran  and  C through  the  interface  proposed 
in  the  Sparse  BLAS  Toolkit[l].  Also  see  the  companion  paper  [2]. 

In  addition  to  the  Sparse  BLAS  Toolkit  interface,  developers  have  access  to  lightweight  kernel  routines. 
These  Sparse  BLAS  Lite  routines  are  unique  to  each  parameter  combination  of  the  higher-level  Toolkit 
interface.  The  Lite  routines  are  designed  for  minimal  overhead;  they  have  no  case  statements,  nor 
elaborate  error-detection  overhead.  Thus,  they  are  ideal  for  use  on  small  matrices  or  to  be  used  as 
efficient  building  blocks  in  higher-level  routines.  Some  typical  examples  of  the  Lite  routines: 


C <-  A’  * B 
C <-  A * B + C 
C <-  alpha*A*B  + b*C 
C <-  D*A~ (-1) *B  + C 


CSR_MatMult_CATB_double() 
CSR_MatMult_CABC_double() 
CSR_MatMult _CaBbC_double ( ) 
CSR_MatTriangSlvLD_CDABC_double ( ) 


These  lightweight  kernel  routines  are  generated  from  a small  number  of  source  lines  (less  than  5000  for 
the  storage  formats  currently  supported)  by  defining  and  expanding  macros  for  successively  restrictive 
sets  of  calling  sequence  parameters.  This  allows  changes  to  the  core  source  code,  made  for  optimization 
or  debugging,  to  be  rapidly  and  automatically  propagated  to  all  affected  kernel  routines  (approximately 
130,000  lines  of  code). 

Section  2 gives  an  introduction  to  the  source  code  generation  mechanism  for  Toolkit’s  underlying  “Lite” 
kernel  library.  Section  3 provides  interface  specifications  for  the  Toolkit  routines  provided  in  this  release, 
and  Section  4 gives  the  function  prototypes  for  the  “Lite”  interface  routines  for  the  VBR  (variable  block 
row)  format  as  an  example.  Prototypes  for  other  formats  are  similar,  and  can  be  obtained  directly 
from  the  header  files  in  the  include  subdirectory.  Installation  instructions  for  the  library  are  provided 
in  Appendix  A. 
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Note  for  Fortran  Users:  The  interfaces  described  in  this  user  guide  are  C interfaces.  For  the  Toolkit 
layer,  the  Fortran  interfaces  are  similar,  except  that  all  arguments  are  passed  by  reference  (that  is, 
typical  Fortran  style).  The  “Lite”  interface  is  not  currently  available  in  Fortran,  because  of  the  need  for 
long  routine  names.  This  restriction  will  be  re-examined  for  future  releases. 


2 Source  Code  Generation 

The  SRC -GEN  directory  contains  the  following  generic  source  code  files. 


bcomm. c 
bscmm . c 
bscmts . c 


bsrmm . c 
bsrmts . c 
coomm . c 


cscnun.  c 
cscmts . c 
csrmm . c 


csrmts . c 
vbrmm . c 
vbrmts . c 


Also  provided  are  generator  scripts  for  creating  the  NIST  Sparse  BLAS  kernel  routines  from  these  generic 
source  files. 

These  source  files  are  used  as  “master  files",  and  are  written  in  such  a way  that  special  case  routines 
can  be  generated  by  relatively  simple  shell  scripts  which  use  “sed"  and  “awk"  for  text  replacement. 
The  approach  saves  considerable  programming  effort  by  generating  most  source  files  automatically,  and 
reduces  errors  by  insuring  that  any  changes  are  propagated  throughout  all  of  the  related  source  code. 

The  master  files  provide  working  source  code  for  the  most  general  version  of  the  kernel  routine.  This  is 
where  real  programming  effort  should  be  expended  to  optimized  the  library.  The  code  is  commented  with 
tags  which  can  be  used  to  selectively  delete  code  for  special  case  routines.  The  "rules"  for  creating  each 
special  case  file  are  defined  in  the  SRC_GEN/kernels  subdirectory.  The  kernels  subdirectory  contains 
the  files 


CAB 

CADBbC 

CDADBC 

CaADB 

CaDABbC 

CABC 

CDAB 

CDADBbC 

CaADBC 

CaDADB 

CABbC 

CDABC 

CaAB 

CaADBbC 

CaDADBC 

CADB 

CDABbC 

CaABC 

CaDAB 

CaDADBbC 

CADBC 

CDADB 

CaABbC 

CaDABC 

one  representing  each  of  the  specializations  from  the  generic  master  code,  along  with  kernel  files  for 
the  master  codes.  Each  of  these  kernel  files  contains  pointers  to  appropriate  "Definition"  files,  in  the 
directory  SRC-GEN/Defs,  which  are  used  to  build  up  the  sed  script  for  the  text  replacement  to  generate 
the  kernel  routines. 

For  typical  use,  these  kernel  and  definition  files  would  never  have  to  be  touched.  Many  modifications 
(say  for  optimization)  can  be  made  to  the  master  source  files  without  requiring  any  change  whatsoever 
to  the  file  generation  mechanism.  The  only  source  code  changes  which  would  affect  code  generation 
would  be  those  which  alter  the  relationship  between  the  comment  tags  and  the  related  source.  A more 
detailed  explanation  of  the  mechanism,  and  requirements  for  modifications,  will  be  forthcoming  in  the 
1.0  release. 
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NIST  Sparse  BLAS 


After  making  any  necessary  changes  to  these  "master"  source  files,  the  library  source  files  may  be 
generated  via  the  "create"  script  (automated  in  the  "make"  process  in  this  directory  with  make  install 
or  "make  re-install"). 

**  IMPORTANT  NOTE  ** 

Any  changes  to  source  for  any  routines  below  the  Toolkit  interface  layer  MUST  be  made  in  the  ../S- 
RC-GEN  directory  to  be  retained  and  propagated  to  all  appropriate  kernel  routines.  Changes  to  the 
Toolkit  interface  routines,  however,  should  be  made  directly  in  the  directory  . . /src_tk  [c  1 f ] .) 
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Toolkit  Interface  Descriptions 


Version  0.9 


4 


July  1996 


dbcomm 


Sparse  BLAS  Toolkit  Interface 


dbcomm 


Name  dbcomm 

Calling  Sequence 

void  dbcommC  const  int  transa,  const  int  mb,  const  int  n,  const  int  kb, 
const  double  alpha,  const  int  descra[] , const  double  val [] , 
const  int  bindx[],  const  int  bjndx[],  const  int  bnnz, 
const  int  lb,  const  double  b[],  const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc, 
double  work [] , const  int  lwork) ; 


Functionality 


Block  coordinate  format  matrix-matrix  multiply. 


C <-  qAB  + fiC 
C <r-  aAr  B + 0C 


Arguments 


int  transa 

Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 

int  mb 

Number  of  block  rows  in  matrix  A 

int  n 

Number  of  columns  in  matrix  c 

int  kb 

Number  of  block  columns  in  matrix  A 

double  alpha 

Scalar  parameter 

double  beta 

Scalar  parameter 

int  descra[] 

Descriptor  argument.  Five  element  integer  array 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 
0 : unknown 


0.9 
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dbcomm 


Sparse  BLAS  Toolkit  Interface 


dbcomm 


double  *val 
int  *bindx 

int  *bjndx 

int  bnnz 
int  lb 
double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 


1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  bnnz  consisting  of  the  block  row 
indices  of  the  entries  of  A. 

integer  array  of  length  bnnz  consisting  of  the  block  column 
indices  of  the  entries  of  A. 
number  of  block  entries 
dimension  of  blocks 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 


0.9 
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dbscmm 


Sparse  BLAS  Toolkit  Interface 


dbscmm 


Name  dbscmm 

Calling  Sequence 


void  dbscmm ( const  int  transa,  const  int  mb,  const  int  n,  const  int  kb, 
const  double  alpha,  const  int  descra[]  , const  double  val[], 
const  int  bindx[],  const  int  bpntrb[],  const  int  bpntre[], 
const  int  lb,  const  double  b[],  const  int  ldb, 
const  double  beta,  double  c [] , const  int  ldc, 
double  work  [] , const  int  Iwork) ; 


Functionality 


Block  sparse  column  format  matrix-matrix  multiply. 


C <-  aAB  + 0C 
C f-  ccAr B + 0C 


Arguments 


int  transa 


int  mb 
int  n 
int  kb 

double  alpha 
double  beta 
int  descraQ 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  block  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Number  of  block  columns  in  matrix  A 
Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descrafO]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descrafl]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 
0 : unknown 
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dbscmm 


Sparse  BLAS  Toolkit  Interface 


dbscmm 


double  *val 
int  *bindx 

int  *bpntrb 

int  *bpntre 

int  lb 
double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 


1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  bnnz  consisting  of  the  block  row 
indices  of  the  entries  of  A. 

integer  array  of  length  mb  such  that  bpntrb(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  first  block  entry  of 
the  j-th  column  of  A. 

integer  array  of  length  mb  such  that  bpntre(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  last  block  entry  of 
the  j-th  column  of  A. 
dimension  of  blocks 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 


0.9 
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dbscsm 


Sparse  BLAS  Toolkit  Interface 


dbscsm 


Name  dbscsm 

Calling  Sequence 


void  dbscsm(  const  int  transa,  const  int  mb,  const  int  n, 
const  int  unitd,  const  double  dv[]  , 

const  double  alpha,  const  int  descraf] , const  double  val [] , 
const  int  bindx[],  const  int  bpntrb[],  const  int  bpntre[], 
const  int  lb,  const  double  b[] , const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc , 
double  work[]  , const  int  lwork) ; 


Functionality 


Block  sparse  column  format  triangular  solve. 


C <-  aDA~l  B + PC  C «-  aDA~TB  + pC 
C aA-'DB  + PC  C <-  aA~TDB  + pC 


Arguments 


int  transa 


int  mb 
int  n 
int  unitd 


double  alpha 
double  beta 
int  descraf] 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  block  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Type  of  scaling: 

1 : Identity  matrix  (argument  dv[]  is  ignored) 

2 : Scale  on  left  (row  scaling) 

3 : Scale  on  right  (column  scaling) 

Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew ( Anti )-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 


0.9 
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dbscsm 


Sparse  BLAS  Toolkit  Interface 


dbscsm 


double  *val 
int  *bindx 

int  *bpntrb 

int  *bpntre 

int  bnnz 
int  lb 
double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 


0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 

0 : unknown 

1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  bnnz  consisting  of  the  block  row 
indices  of  the  entries  of  A. 

integer  array  of  length  mb  such  that  bpnt,rb(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  first  block  entry  of 
the  j-th  column  of  A. 

integer  array  of  length  mb  such  that  bpntre(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  last  block  entry  of 
the  j-th  column  of  A. 
number  of  block  entries 
dimension  of  blocks 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 


0.9 
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dbsrmm 


Sparse  BLAS  Toolkit  Interface 


dbsrmm 


Name  dbsrmm 

Calling  Sequence 


void  dbsrmm(  const  int  transa,  const  int  mb,  const  int  n,  const  int  kb, 
const  double  alpha,  const  int  descra[],  const  double  val [] , 
const  int  bindx[],  const  int  bpntrb  []  , const  int  bpntre[], 
const  int  lb,  const  double  b[],  const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc, 
double  work  []  , const  int  lwork) ; 


Functionality 


Block  sparse  row  format  matrix-matrix  multiply. 


C <-  aAB  + pC 
C +-  aArB  + 0C 


Arguments 


int  transa 


int  mb 
int  n 
int  kb 

double  alpha 
double  beta 
int  descraf] 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  block  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Number  of  block  columns  in  matrix  A 
Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 
0 : unknown 


0.9 
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dbsrmm 


Sparse  BLAS  Toolkit  Interface 


dbsrmm 


double  *val 
int  * bindx 

int  *bpntrb 

int  *bpntre 

int  lb 
double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 


1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  bnnz  consisting  of  the  block  column 
indices  of  the  entries  of  A. 

integer  array  of  length  mb  such  that  bpntrb(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  first  block  entry  of 
the  j-th  row  of  A. 

integer  array  of  length  mb  such  that  bpntre(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  last  block  entry  of 
the  j-th  row  of  A. 
dimension  of  blocks 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 


0.9 
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dbsrsm 


Sparse  BLAS  Toolkit  Interface 


dbsrsm 


Name  dbsrsm 

Calling  Sequence 


void  dbsrsm(  const  int  transa,  const  int  mb,  const  int  n, 
const  int  unitd,  const  double  dv [] , 

const  double  alpha,  const  int  descra[],  const  double  val [] , 
const  int  bindx[],  const  int  bpntrb[],  const  int  bpntre[], 
const  int  lb,  const  double  b[] , const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc , 
double  work[] , const  int  lwork) ; 


Functionality 


Block  sparse  row  format  triangular  solve. 

C *r-  otDA~l B + pC  C <-  otDA-rB  + pC 

C «—  aA~x  DB  + (3C  C <—  qA~t  DB  + 0C 


Arguments  int  transa 

int  mb 
int  n 
int  unitd 


double  alpha 
double  beta 
int  descra[] 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  block  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Type  of  scaling: 

1 : Identity  matrix  (argument  dv[]  is  ignored) 

2 : Scale  on  left  (row  scaling) 

3 : Scale  on  right  (column  scaling) 

Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew ( Anti )-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 


0.9 
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dbsrsm 


Sparse  BLAS  Toolkit  Interface 


dbsrsm 


double  *val 
int  *bindx 

int  *bpntrb 


int  *bpntre 


int  lb 
double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 


0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 

0 : unknown 

1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  bnnz  consisting  of  the  block  column 
indices  of  the  entries  of  A. 

integer  array  of  length  mb  such  that  bpntrb(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  first  block  entry  of 
the  j-th  row  of  A. 

integer  array  of  length  mb  such  that  bpntre(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  last  block  entry  of 
the  j-th  row  of  A. 
dimension  of  blocks 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 


0.9 


14 


dcoomm 


Sparse  BLAS  Toolkit  Interface 


dcoomm 


Name  dcoomm 

Calling  Sequence 

void  dcoomm ( const  int  transa,  const  int  m,  const  int  n,  const  int  k, 


Functionality 

const  double  alpha,  const  int  descraG  , const  double  val [] 
const  int  indx[],  const  int  jndx[],  const  int  nnz, 
const  double  b[] , const  int  ldb, 
const  double  beta,  double  c [] , const  int  ldc, 
double  work[],  const  int  lwork) ; 

Coordinate  format  matrix-matrix  multiply. 


C +-  aAB  + PC 

C <-  aAT B + PC 

Arguments  int  transa 

Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

int  m 

int  n 

int  k 

double  alpha 
double  beta 
int  descra[] 

1 : operate  with  transpose  matrix 

Number  of  rows  in  matrix  A 

Number  of  columns  in  matrix  c 

Number  of  columns  in  matrix  A 

Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 

0 : unknown 
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dcoomm 


dcoomm 


Sparse  BLAS  Toolkit  Interface 


double  *val 
int  *indx 
int  *jndx 
int  nnz 
double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 


1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  nnz  containing  row  indices 
integer  array  of  length  nnz  containing  column  indices 
Number  of  nonzero  matrix  entries 
rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
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dcscmm 


Sparse  BLAS  Toolkit  Interface 


dcscmm 


Name  dcscmm 

Calling  Sequence 

void  dcscmm ( const  int  transa,  const  int  m,  const  int  n,  const  int  k, 

const  double  alpha,  const  int  descra[]  , const  double  val[], 
const  int  indx[],  const  int  pntrb[],  const  int  pntre[], 
const  double  b[] , const  int  ldb, 
const  double  beta,  double  c [] , const  int  ldc, 
double  work[] , const  int  lwork) ; 


Functionality 


Compressed  sparse  column  format  matrix-matrix  multiply. 


C +-  OiAB  + (3C 
C <-  ocAtB  + 0C 


Arguments 


descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 
0 : unknown 


int  transa 


int  m 
int  n 
int  k 

double  alpha 
double  beta 
int  descra[] 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Number  of  columns  in  matrix  A 
Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
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dcscmm 


Sparse  BLAS  Toolkit  Interface 


dcscmm 


double  *val 
int  *indx 
int  *pntrb 

int  *pntre 

double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 


1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  nnz  containing  row  indices 
integer  array  of  length  k such  that  pntrb(j)-pntrb(l) 
points  to  location  in  val  of  the  first  nonzero  element 
in  column  j 

integer  array  of  length  k such  that  pntre(j)-pntrb(  1 ) 
points  to  location  in  val  of  the  last  nonzero  element 
in  column  j 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lw'ork.  lwork  should  be  at  least  max(m,n) 
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dcscsm 


Sparse  BLAS  Toolkit  Interface 


dcscsm 


Name  dcscsm 

Calling  Sequence 


void  dcscsm(  const  int  transa,  const  int  m,  const  int  n, 
const  int  unitd,  const  double  dv  []  , 

const  double  alpha,  const  int  descra[] , const  double  val  [] , 
const  int  indx[],  const  int  pntrb[],  const  int  pntre[], 
const  double  b[] , const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc, 
double  work  [] , const  int  lwork) ; 


Functionality 


Compressed  sparse  column  format  triangular  solve. 


C <r-  aDA~lB  + pC  C <-  aDA~TB  + PC 

C <-  aA~l  DB  + PC  C <r-  aA~TDB  + pC 

Arguments 


descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 


int  transa  Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 

int  m Number  of  rows  in  matrix  A 

int  n Number  of  columns  in  matrix  c 

int  unitd  Type  of  scaling: 

1 : Identity  matrix  (argument  dv|]  is  ignored) 

2 : Scale  on  left  (row  scaling) 

3 : Scale  on  right  (column  scaling) 

double  alpha  Scalar  parameter 

double  beta  Scalar  parameter 

int  descraQ  Descriptor  argument.  Five  element  integer  array: 
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dcscsm 


Sparse  BLAS  Toolkit  Interface 


dcscsm 


double  *val 
int  *indx 
int  *pntrb 


int  *pntre 


double  *b 
int  ldb 
double  *c 
int  Idc 

double  *work 
int  lwork 


0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 

0 : unknown 

1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  nnz  containing  row  indices 
integer  array  of  length  k such  that  pntrb(j)-pntrb(l) 
points  to  location  in  val  of  the  first  nonzero  element 
in  column  j 

integer  array  of  length  k such  that  pntre(j)-pnt.rb(l) 
points  to  location  in  val  of  the  last  nonzero  element 
in  column  j 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 
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dcsrmm 


Sparse  BLAS  Toolkit  Interface 


dcsrmm 


Name  dcsrmm 

Calling  Sequence 

void  dcsrmm(  const  int  transa,  const  int  m,  const  int  n,  const  int  k, 

const  double  alpha,  const  int  descraf] , const  double  val [] , 
const  int  indx[],  const  int  pntrb[],  const  int  pntre  []  , 
const  double  b[],  const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc, 
double  work  [] , const  int  lwork) ; 


Functionality 


Compressed  sparse  row  format  matrix-matrix  multiply. 


C +-  aAB  + pC 
C <-  aATB  + PC 


Arguments 


descrafO]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 
0 : unknown 


int  transa 


int  m 
int  n 
int  k 

double  alpha 
double  beta 
int  descraf] 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Number  of  columns  in  matrix  A 
Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
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dcsrmm 


Sparse  BLAS  Toolkit  Interface 


dcsrmm 


double  *val 
int  *indx 
int  *pntrb 

int  *pntre 

double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 


1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  nnz  containing  column  indices 
integer  array  of  length  k such  that  pntrb(j)-pntrb(l) 
points  to  location  in  val  of  the  first  nonzero  element 
in  row  j 

integer  array  of  length  k such  that  pntre(j)-pntrb(l) 
points  to  location  in  val  of  the  last  nonzero  element 
in  row  j 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
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dcsrsm 


Sparse  BLAS  Toolkit  Interface 


dcsrsm 


Name  dcsrsm 

Calling  Sequence 

void  dcsrsmC  const  int  transa,  const  int  m,  const  int  n 


Functionality 

const  int  unitd,  const  double  dv[], 

const  double  alpha,  const  int  descra[] , const  double  val [] 
const  int  indx[],  const  int  pntrb[]  , const  int  pntre[], 
const  double  b[],  const  int  ldb, 
const  double  beta,  double  c[],  const  int  Idc, 
double  work[] , const  int  lwork) ; 

Compressed  sparse  row  format  triangular  solve. 


C <-  aDA~lB  + pC  C aDA-TB  + PC 

C <-  a.A~x  DB  + PC  C <-  aA~TDB  + pC 

Arguments  int  transa 

Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 

int  m 

int  n 

int  unitd 

Number  of  rows  in  matrix  A 

Number  of  columns  in  matrix  c 

Type  of  scaling: 

double  alpha 
double  beta 
int  descraQ 

1 : Identity  matrix  (argument  dv[]  is  ignored) 

2 : Scale  on  left  (row  scaling) 

3 : Scale  on  right  (column  scaling) 

Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 
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dcsrsm 


Sparse  BLAS  Toolkit  Interface 


dcsrsm 


double  *val 
int  *indx 
int  *pntrb 


int  *pntre 


double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 


0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 

0 : unknown 

1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  nnz  containing  column  indices 
integer  array  of  length  k such  that  pntrb(j)-pntrb(l) 
points  to  location  in  val  of  the  first  nonzero  element 
in  row  j 

integer  array  of  length  k such  that  pntre(j)-pntrb(l) 
points  to  location  in  val  of  the  last  nonzero  element 
in  row  j 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  max(m,n) 
length  of  work  array 
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dvbrmm 


Sparse  BLAS  Toolkit  Interface 


dvbrmm 


Name  dvbrmm 

Calling  Sequence 


void  dvbrmm ( const  int  transa,  const  int  mb,  const  int  n,  const  int  kb, 
const  double  alpha,  const  int  descra[] , const  double  val [] , 
const  int  indx[],  const  int  bindx[], 
const  int  rpntr[],  const  int  cpntr[], 
const  int  bpntrbf]  , const  int  bpntre[], 
const  double  b[],  const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc , 
double  work  [] , const  int  lwork) ; 


Functionality 


Variable  block  row  format  matrix-matrix  multiply. 


C <-  aAB  + 0C 
C +-  aATB  + (3C 


Arguments 


int  transa 


int  mb 
int  n 
int  kb 

double  alpha 
double  beta 
int  descraf] 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  block  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Number  of  block  columns  in  matrix  A 
Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descra[l]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 
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dvbrmm 


Sparse  BLAS  Toolkit  Interface 


dvbrmm 


s 

descra[4]  repeated  indice  (not  currently  supported) 

0 : unknown 

double  *val 
int  *indx 

1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 
integer  array  of  length  bnnz+1  such  that  the  i-th 
element  of  indx[]  points  to  the  location  in  val  of 
the  (1,1)  element  of  the  i-th  block  entry. 

int,  * bindx 

integer  array  of  length  bnnz  consisting  of  the  block  column 
indices  of  the  entries  of  A. 

int  *rpntr 

integer  array  of  length  mb+1  such  that  rpntr(i)-rpntr(l) 
is  the  row  index  of  the  first  point  row  in  the  i-th 
block  row.  rpntr(mb-fl)  is  set  to  m+rpntr(l). 

Thus,  the  number  of  point  rows  in  the  i-th  block  row  is 

int  *cpntr 

rpntr(i+l)-rpntr(i). 

integer  array  of  length  kb+1  such  that  cpntr(j)-cpntr(l) 
is  the  column  index  of  the  first  point  column  in  the  j-th 
block  column,  cpntr(kb-t-l)  is  set  to  k+cpntr(l). 

Thus,  the  number  of  point  columns  in  the  j-th  block  column 
is  cpntr(j  + l)-cpntr(j). 

int  *bpntrb 

integer  array  of  length  mb  such  that  bpntrb(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  first  block  entry  of 
the  j-th  row  of  A. 

int  *bpntre 

integer  array  of  length  mb  such  that  bpntre(i)-bpntrb(l) 
points  to  location  in  bindx  of  the  last  block  entry  of 
the  j-th  row  of  A. 

double  *b 
int  ldb 
double  *c 
int  ldc 

double  *work 
int  lwork 

rectangular  array  with  leading  dimension  ldb 
leading  dimension  of  b 

rectangular  array  with  leading  dimension  ldc 
leading  dimension  of  c 

scratch  array  of  length  lwork.  lwork  should  be  at  least  m*n  + max(blocksize)2 
length  of  work  array 
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dvbrsm 


Sparse  BLAS  Toolkit  Interface 


dvbrsm 


Name  dvbrsm 

Calling  Sequence 

void  dvbrsm(  const  int  transa,  const  int  mb,  const  int  n, 
const  int  unitd,  const  double  dv[] , 

const  double  alpha,  const  int  descra[] , const  double  val [] , 
const  int  indx[],  const  int  bindx[],  const  int  rpntr[], 
const  int  cpntr[],  const  int  bpntrb[],  const  int  bpntre[], 
const  double  b[],  const  int  ldb, 
const  double  beta,  double  c[],  const  int  ldc, 
double  work[] , const  int  lwork) ; 


Functionality 


Variable  block  row  format  triangular  solve. 


C <—  aDA~l B + PC  C <-  otDA~r B + (3C 

C 4-  aA~1  DB  + 0C  C aA~T DB  + pC 


Arguments 


int  transa 

int  mb 
int  n 
int  unitd 


double  alpha 
double  beta 
int  descraQ 


Indicates  how  to  operate  with  the  sparse  matrix 

0 : operate  with  matrix 

1 : operate  with  transpose  matrix 
Number  of  block  rows  in  matrix  A 
Number  of  columns  in  matrix  c 
Type  of  scaling: 

1 : Identity  matrix  (argument  dv[]  is  ignored) 

2 : Scale  on  left  (row  scaling) 

3 : Scale  on  right  (column  scaling) 

Scalar  parameter 

Scalar  parameter 

Descriptor  argument.  Five  element  integer  array: 
descra[0]  matrix  structure 

0 : general 

1 : symmetric 

2 : Hermitian 

3 : Triangular 

4 : Skew(Anti)-Symmetric 

5 : Diagonal 

descrajl]  upper/lower  triangular  indicator 

1 : lower 

2 : upper 

descra[2]  main  diagonal  type 

0 : non-unit 

1 : unit 
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dvbrsm 


Sparse  BLAS  Toolkit  Interface 


dvbrsm 


double  *val 

descra[3]  Array  base 

0 : C/C++  compatible 

1 : Fortran  compatible 

descra[4]  repeated  indices  (not  currently  supported) 

0 : unknown 

1 : no  repeated  indices 

scalar  array  of  length  nnz  containing  matrix  entries 

int  *indx 

integer  array  of  length  bnnz+1  such  that  the  i-th  element 

int  *bindx 

of  indx[]  points  to  the  location  in  val  of  the  (1,1)  element 
of  the  i-th  block  entry. 

integer  array  of  length  bnnz  consisting  of  the  block  column 

int  *rpntr 

indices  of  the  entries  of  A. 

integer  array  of  length  mb+1  such  that  rpntr(i)-rpntr(l) 

int  *cpntr 

is  the  row  index  of  the  first  point  row  in  the  i-th  block 
row.  rpntr(mb+l)  is  set  to  m+rpntr(l). 

Thus,  the  number  of  point  rows  in  the  i-th  block  row  is 
rpntr(i+l)-rpntr(i). 

integer  array  of  length  kb+1  such  that  cpntr(j)-cpntr(l) 

int  *bpntrb 

is  the  column  index  of  the  first  point  column  in  the  j-th 
block  column.  cpntr(kb+l)  is  set  to  k+cpntr(l). 

Thus,  the  number  of  point  columns  in  the  j-th  block  column  is 
cpntr  ( j + 1 ) -cpntr  ( j ) . 

integer  array  of  length  mb  such  that  bpnt.rb(i)-bpntrb(l) 

int  *bpntre 

points  to  location  in  bindx  of  the  first  block  entry  of 
the  j-th  row  of  A. 

integer  array  of  length  mb  such  that  bpntre(i)-bpntrb(l) 

double  *b 

points  to  location  in  bindx  of  the  last  block  entry  of 
the  j-th  row  of  A. 

rectangular  array  with  leading  dimension  ldb 

int  ldb 

leading  dimension  of  b 

double  *c 

rectangular  array  with  leading  dimension  ldc 

int  ldc 

leading  dimension  of  c 

double  *work 

scratch  array  of  length  lwork.  lwork  should  be  at  least  m*n  + 

int  lwork 

length  of  work  array 

max(blocksize)2 


| 
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4 


uLite"  Function  Prototypes:  VBR  Example 


This  section  is  provided  to  give  an  illustration  of  the  naming  convention  and  corresponding  prototypes 
for  the  automatically  generating  lightweight  functions.  Since  there  are  over  250  functions  for  the  VBR 
storage  format  alone,  it  is  important  that  these  functions  and  prototypes  follow  a predictable  scheme. 

If  a user  is  familiar  with  the  prototype  for  the  most  generic  function  (corresponding  to  the  matrix 
matrix  multiplication  kernel  CaABbC  and  the  matrix  triangular  solve  kernel  CaDADBbC),  then  special 
case  prototypes  can  be  predicted  by  changing  the  kernel  name  and  dropping  out  any  corresponding 
unnecessary  arguments  from  the  argument  list.  For  example,  the  vector  version  of  the  routine 

void  VBR_MatTriangSlvLU_CaDADBbC_double(  const  int  mb,  const  int  n,  const  double 
*dvl , const  double  *dvr , const  double  alpha,  const  double  *val , const  int  *indx,  const  int 
*bindx , const  int  *rpntr , const  int  *cpntr , const  int  *bpntrb , const  int  *bpntre , const 
double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind_base) ; 

can  be  obtained  by  changing  the  function  name  to  VBR,  VecTriangSlvLU. CaDADBbC -double, 
and  eliminating  the  arguments  { n,  ldb,  ldc,  to  arrive  at  the  prototype: 

void  VBR_VecTriangSlvLU_CaDADBbC -double ( const  int  mb,  const  double  *dvl , const 
double  *dvr , const  double  alpha,  const  double  *val , const  int  *indx , const  int  *bindx , const 
int  *rpntr , const  int  *cpntr , const  int  *bpntrb,  const  int  *bpntre , const  double  *b , const 
double  beta,  double  *c , double  *work , const  int  ind_base) ; 

The  protoypes  listed  in  this  section  further  illustrate  the  use  of  this  convention. 


4.1  Variable  Block  Row  Matrix  Multiply  Routines 


void  VBR_MatMult_CAB-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  const  int  md-base); 

void  VBR_MatMult_CATB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRsymm_MatMult_CAB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  double  *c,  const  int  ldc,  const  mt  ind-base); 

void  VBRskew _M.atMult_CAB-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int 
*inctx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBRskew -MatMult_CATB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult_CaAB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const  double 
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*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR-MatMult_CaATB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBRsymm_MatMult_CaAB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBRskew _MatMu!t_CaAB_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRskew  JVIat  Mult  _CaATB-doub!e(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult-CABC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult_CATBC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRsymm JMatM«lt_CABC-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRskew_MatMult _CABC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRskew _MatMult_CATBC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult-CaABC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult_CaATBC-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRsymmi-MatMult_CaABC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBRskew-MatMult-CaABC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
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*bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBRskew -MatMult_CaATBC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBR_MatMult_CABbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult_CATBbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBRsymm JVIatMult_CABbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRskew_MatMult_CABbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRskew_MatMult_CATBbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult_CaABbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatMult_CaATBbC-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bmdx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBRsymm_MatMult_CaABbC-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBRskew JMatMult_CaABbC-double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  mt  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBRskew -Mat  Mult_CaATBbC_double(  const  int  mb,  const  int  n,  const  int  kb,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  const  int  ind_base); 


4.2  Variable  Block  Row  Vector  Multiply  Routines 


void  VBR_VecMult_CAB_double(  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const 
int  ind_base); 
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void  VBR_VecMult_CATB_double(  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const 
int  ind.base); 

void  VBRsymm.VecMult.CAB.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const 
int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
const  int  ind.base); 

void  VBRskew_VecMult_CAB .double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const 
int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
const  int  ind.base); 

void  VBRskew.VecMult.CATB  .double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double 
*c,  const  int  ind.base); 

void  VBR.VecMult.CaAB. double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  const  int  ind.base); 

void  VBR.VecMult.CaATB .double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  const  int  ind.base); 

void  VBRsymm.VecMult.CaAB.double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  double  *c,  const  int  ind.base); 

void  VBRskew.VecMult.CaAB. double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  double  *c,  const  int  ind_base); 

void  VBRskew.VecMult.CaATB.double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  double  *c,  const  int  ind.base); 

void  VBR.VecMult.C  ABC  .double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const 
int  ind.base); 

void  VBR.VecMult.C ATBC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const 
int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
const  int  ind.base); 

void  VBRsymm.VecMult.CABC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double 
*c,  const  int  ind.base); 

void  VBRskew.VecMult.CABC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double 
*c,  const  int  ind.base); 


void  VBRskew.VecMult.CATBC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
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const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double 
*c,  const  int  ind_base); 

void  VBR_VecMult_CaABC_double(  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  double  *c,  const  int  ind-base); 

void  VBR.VecMult.CaATBC.double)  const  int  mb,  const  int,  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  double  *c,  const  int  ind.base); 

void  VBRsymm_VecMult_CaABC_double(  const  int  mb,  const  int  kb,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  double  *c,  const  int  ind.base); 

void  VBRskew_VecMult_CaABC_double(  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  double  *c,  const  int  ind.base); 

void  VBRskew_VecMult_CaATBC_double(  const  int  mb,  const  int  kb,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  double  *c,  const  int  ind.base); 

void  VBR.VecMult.CABbC .double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const 
int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double 
beta,  double  *c,  const  int  ind.base); 

void  VBR.VecMult.CATBbC .double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx,  const 
int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double 
beta,  double  *c,  const  int  ind.base); 

void  VBRsymm.VecMult.CABbC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
double  beta,  double  *c,  const  int  ind.base); 

void  VBRskew.VecMult.CABbC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
double  beta,  double  *c,  const  int  ind.base); 

void  VBRskew.VecMult.CATBbC.double)  const  int  mb,  const  int  kb,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
double  beta,  double  *c,  const  int  ind.base); 

void  VBR.VecMult.CaABbC.double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  const  double  beta,  double  *c,  const  int  ind.base); 

void  VBR.VecMult.CaATBbC.double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double  *val, 
const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const 
double  *b,  const  double  beta,  double  *c,  const  int  ind.base); 

void  VBRsymm.VecMult.CaABbC.double)  const  int  mb,  const  int  kb,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
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const  double  *b,  const  double  beta,  double  *c,  const  int  ind-base); 

void  VBRskew.VecMult_CaABbC_double(  const  int  mb,  const  int  kb,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  double  beta,  double  *c,  const  int  ind-base); 

void  VBRskew_VecMult_CaATBbC_doubie(  const  int  mb,  const  int  kb,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  double  beta,  double  *c,  const  int  ind.base); 


4.3  Variable  Block  Row  Matrix  Triangular  Solve  Routines 

void  VBR_MatTriangSllvLU_CAB_doubIe(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CAB_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBR_MatTriangSlvLU_CATB_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatTriangSlvUU-CATB-double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR-MatTriangSlvLU_CaAB_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaAB_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatTriamgSlvLU-CaATB-do«ble(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaATB_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind_base); 

void  VBR_MatTriangSlvLU_CABC-double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CABC_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 
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void  VBR_MatTriangSlvLU_CATBC_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CATBC_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CaABC_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaABC_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CaATBC-double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaATBC_doub!e(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR-MatTriangSlvLU_CABbC_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CABbC_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int  *indx, 
const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const 
int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind_base); 

void  VBR_MatTriangSlvLU_CATBbC-double(  const  int  mb,  const  int  n,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CATBbC_double(  const  int  mb,  const  int  n,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU-CaABbC-double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUB_CaABbC_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU-CaATBbC_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaATBbC_double(  const  int  mb,  const  int  n,  const  double  alpha,  const  double 


35 


*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  Idb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CDAB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int,  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR _MatTriangSlvUU_CBAB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind_base); 

void  VBR_MatTriangSlvLU-CDATB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CDATB-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CaOAB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUB_CaBAB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR-MatTriangSlvLU_CaDATB.double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaDATB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CBABC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind_base); 

void  VBR_MatTriangSlvUU_CDABC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  +bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CDATBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUB_CDATBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU-CaDABC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
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const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaDABC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CaDATBC_double(  const  int  mb,  const  int  n,  const  double  *dvi,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaDATBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CDABbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CDABbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CBATBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU-CDATBbC-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangS!vLU_CaDABbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bmdx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind-base); 

void  VBR_MatTriangSlvUU_CaDABbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind_base); 

void  VBR_MatTriamgSlvLU_CaDATBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaDATBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CADB_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind-base); 
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void  VBR_MatTriangSlvUU_CADB_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CATBB_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBR_MatTriangSlvUB_CATBB_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBRJVIatTriangSlvLU.Ca  ABB  .double)  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBRJVIatTriangSlvBU_CaABB_doMb!e(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBR_MatTriaogSlvLB_CaATBB_doubie(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBR_MatTriangSlvUU_CaATBB_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  const  int  ind.base); 

void  VBR-MatTrIamgSlvLU_CABBC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvBU_CAOBC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CATBBC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangS!vUU_CATDBC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBRJVIatTriangSlvLU.CaADBC.double)  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvUU_CaABBC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBRJVIatTriangSlvLU.CaATBBC.double)  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
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alpha,  const  double  *val,  const  int  *indx,  const  int  *bmdx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvUU.CaATDBC.double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBRJVIatTriangSlvLU_CADBbC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriamgSlvUU_CADBbC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CATDBbC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvUU_CATDBbC.double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
*val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre, 
const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CaADBbC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind.base); 

void  VBR_MatTriangSIvUU_CaAOBbC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const  double 
alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind.base); 

void  VBR_MatTriangSlvLU.CaATDBbC_double(  const  int  mb,  const  int  n,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind.base); 

void  VBR_MatTriangSlvUB_CaATOBbC_doub!e(  const  int  mb,  const  int  n,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  ’•'cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind.base); 

void  VBR_MatTriangSivLU_CDABB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvUU_CDADB_doub!e(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CDATDB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  ^cpntr,  const  int  *bpntrb, 
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const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvUU_CBATBB-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBRJVtatTriangSlvLU-CaOADB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr, 
const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind.base); 

void  VBR-MatTriangS!vUU_CaDABB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr, 
const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind-base); 

void  VBR_MatTriangSlvLU_CaBATBB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr, 
const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind-base); 

void  VBR_MatTriangSlvUB-CaBATBB_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr, 
const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind-base); 

void  VBR-MatTriamgSlvLU_CDADBC-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CBABBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvLU_CBATBBC-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  ’"bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvUU_CDATDBC-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const  int  ind.base); 

void  VBR_MatTriangSlvLU_CaBABBC.double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr, 
const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind.base); 

void  VBR_MatTriangSlvUU_CaJ3ABBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr, 
const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind.base); 

void  VBR_MatTriangSlvLU_CaBATBBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
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double  *dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const 
int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaDATOBC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const 
int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind_base); 

void  VBRJVIatTriangSlvLU_CDADBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind_base); 

void  VBR_MatTriangSIvUU_CDADBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const  double 
*dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb, 
const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double  *work,  const 
int  ind-base); 

void  VBR_MatTriangSlvLU-CDATDBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CDATDBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const  int  ldc,  double 
*work,  const  int  ind_base); 

void  VBR_MatTriangSivLB_CaDABBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bincbc,  const  int  *rpntr,  const  int 
*cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const 
int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangSlvUU_CaDABBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int 
*cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const 
int  ldc,  double  *work,  const  int  ind-base); 

void  VBR_MatTriangS!vLU_CaDATDBbC-double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int 
*cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const 
int  ldc,  double  *work,  const  int  ind_base); 

void  VBR_MatTriangS!vUU_CaBATBBbC_double(  const  int  mb,  const  int  n,  const  double  *dvl,  const 
double  *dvr,  const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int 
*cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  int  ldb,  const  double  beta,  double  *c,  const 
int  ldc,  double  *work,  const  int  ind-base); 
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4A  Variable  Block  Row  Vector  Triangular  Solve  Routines 


void  VBR_VecTriangSlvLO_CAB -double)  const  int  mb,  const  double  *val,  const  int  *indx,  const  int  *bindx, 
const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const  int 
ind-base); 

void  VBR-VecTriamgSlvUB_CAB -double)  const  int  mb,  const  double  *val,  const  int  *indx,  const  int  *bindx, 
const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const  int 
ind-base); 

void  VBR-VecTriangSlvLU_CATB -double)  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const 
int  ind-base); 

void  VBR_VecTriangSlvUU_CATB -double)  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  const 
int  ind-base); 

void  VBR_VecTriamgSlvLU_CaAB_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvUU-CaAB_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CaATB -double)  const  int  mb,  const  double  alpha,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CaATB -double)  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  const  int  ind-base); 

void  VBR-VecTriangSlvLU_CABC_double(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
double  *work,  const  int  ind-base); 

void  VBR-VecTriangSlvUU-CABC-double)  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CATBC_doub!e(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
double  *work,  const  int  ind-base); 

void  VBR_VecTriamgSlvUXJ_CATBC_double(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c, 
double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CaABC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  ’''cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind_base); 
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void  VBR_VecTriangSlvUU_CaABC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU-CaATBC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CaATBC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CABbC-double(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double 
beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvUU_CABbC-double(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double 
beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CATBbC_double(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double 
beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CATBbC_double(  const  int  mb,  const  double  *val,  const  int  *indx,  const  int 
*bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double 
beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CaABbC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CaABbC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CaATBbC-double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CaATBbC_double(  const  int  mb,  const  double  alpha,  const  double  *val,  const 
int  *incbc,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CDAB_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CDAB_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CDATB_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const  int 
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*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvUU.CDATB  .double)  const  int  mb,  const  double  *dvl,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind_base); 

void  VBR.VecTriamgSlvLU.CaOAB  .double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR.VecTriamgSlvBU.CaOAB  .double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangS!vLU_CaDATB_double(  const  int  mb,  const  double  *dvl,  const  double  aJpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriamgSlvBU.CaBATB.double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR.VecTriangSlvLU.CBABC.double)  const  int  mb,  const  double  *dvl,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvBU_CBABC_dotible(  const  int  mb,  const  double  *dvl,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind.base); 

void  VBR-VecTriangSlvLU-CDATBC_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CDATBC_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvLU.CaDABC .double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvUB.CaBABC.double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  ^rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvLU.CaDATBC.double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvUU.CaBATBC.double)  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
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*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CDABbC_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvUU_CDABbC_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CDATBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind_ba.se); 

void  VBR_VecTriangSlvUU_CDATBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CaDABbC_double(  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CaDABbC_doubIe(  const  int  mb,  const  double  *dvl,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CaDATBbC_double(  const  int  mb,  const  double  *dvl,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR-VecTriangSlvUU_CaDATBbC -double)  const  int  mb,  const  double  *dvl,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CADB_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CADB-double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CATDB_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CATDB_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvLU_CaADB_double(  const  int  mb,  const  double  *dvT,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  const  int  ind-base); 
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void  VBR_VecTriangSlvUU_CaADB_double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  const  int  ind_base); 

void  VBR.VecTriangSlvLU_CaATOB_double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  const  int  ind-base); 

void  VBR_VecTriangSlvUU_CaATDB_double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  const  int  ind.base); 

void  VBR.VecTriangSlvLU_CAOBC_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const  int 
*indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double  *b, 
double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU.CABBC_doubIe(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU_CATBBC.double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR.VecTriangSlvUU_CATBBC.doiible(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR-VecTriangSlvLU_CaABBC„double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CaADBC_double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
’"bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU_CaATDBC_double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUB_CaATBBC_double(  const  int  mb,  const  double  *dvr,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvLU_CADBbC_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CADBbC_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLB_CATBBbC_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
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int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvUU_CATDBbC_double(  const  int  mb,  const  double  *dvr,  const  double  *val,  const 
int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int  *bpntre,  const  double 
*b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU.CaADBbC_double(  const  int  mb,  const  double  *dvr,  const  double  alpha,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CaADBbC.double(  const  int  mb,  const  double  *dvr,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU-CaATDBbC.double(  const  int  mb,  const  double  *dvr,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CaATBBbC_double(  const  int  mb,  const  double  *dvr,  const  double  alpha, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU_CDADB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvUU_CDADB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU_CDATDB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bincbc,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriaegSlvUU_CDATDB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU_CaDADB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CaDADB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvLU_CaDATDB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  alpha,  const  double  +vaJ,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU.CaDATDB_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
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*bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CDADBC_doub!e(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CBADBC-double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriangSlvLU_CDATBBC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR_VecTriangSlvUB.CBATBBC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind_base); 

void  VBR_VecTriaegSlvLU_CaBABBC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind-base); 

void  VBR.VecTriamgS!vBB_CaBABBC.double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int 
*bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangSlvLU_CaBATBBC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const 
int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR-VecTriangSlvUU.CaBATBBC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const 
int  *bpntrb,  const  int  *bpntre,  const  double  *b,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR.VecTriangS!vLB_CBABBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU.CDADBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr,  const 
double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const  int 
*bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR-VecTriangSlvLU-CBATBBbC.dooble(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  +bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CBATBBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const  int  *bpntrb,  const 
int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvLU_CaOAOBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const 
int  bpntrb,  const  int  *bpntre,  const  double  +b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 
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void  VBR_VecTriangSlvUU_CaDADBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const 
int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  md.base); 

void  VBR_VecTriangSlvLU_CaDATDBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const 
int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind.base); 

void  VBR_VecTriangSlvUU_CaDATDBbC_double(  const  int  mb,  const  double  *dvl,  const  double  *dvr, 
const  double  alpha,  const  double  *val,  const  int  *indx,  const  int  *bindx,  const  int  *rpntr,  const  int  *cpntr,  const 
int  *bpntrb,  const  int  *bpntre,  const  double  *b,  const  double  beta,  double  *c,  double  *work,  const  int  ind-base); 


A Installation  Instructions 


The  installation  of  the  Sparse  BLAS  Toolkit  is  automated  with  the  “make"  utility.  To  use  “make"  to 
build  the  library: 


1.  Edit  the  file  . /makefile. def  to  reflect  your  system  setup: 

• The  minimum  installation  requires  an  ANSI  C compiler. 

• An  extended  installation  which  includes  Fortran  callable  routines  and  testers  is  available.  If 
the  presence  of  a Fortran  compiler  is  indicated  in  the  makefile. def  file,  the  extended  version 
will  be  installed. 

• The  archival  process  by  default  uses  “ranlib".  If  this  is  not  available  on  your  system,  set 
HASRANLIB  to  T . 

2.  Type: 


make  install  * 

make  installc 
make  library 

make  teste 

make  testf77  * 


to  build  the  library  AND  make  and  run 
the  C and  Fortran  testers 

to  build  the  library  AND  make  and  run  the  C testers 

to  build  the  archive  file  ./lib/libsptk.a 

(tests  are  not  built) 

to  build  and  run  the  C testers 

(library  must  be  pre-built) 

to  build  and  run  the  Fortran  testers 

(library  must  be  pre-built) 


* requires  a Fortran  compiler 


3.  For  space-saving  cleanup,  type  "make  clean"  to  remove  all  .o 
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